ALGEBRAIC PROPERTIES OF A DISORDERED 
ASYMMETRIC GLAUBER MODEL 



ARVIND AYYER 

Abstract. We consider an asymmetric variant of disordered Gla- 
uber dynamics of Ising spins on a one-dimensional lattice, where 
each spin flips according to the relative state of the spin to its 
left. Moreover, each bond allows for two rates; flips which equalize 
nearest neighbor spins, and flips which "unequalize" them. In addi- 
tion, the leftmost spin flips depending on the spin at that site. We 
explicitly calculate all eigenvalues of the transition matrix for all 
system sizes and conjecture a formula for the normalization factor 
of the model. We then analyze two limits of this model, which are 
analogous to ferromagnetic and antiferromagnetic behavior in the 
Ising model for which we are able to prove an analogous formula 
for the normalization factor. 



1. Introduction 

The Ising model with Glauber dynamics [1] has been an extremely 
important source for understanding time dependent behavior in the 
Ising model. Moreover, the simplicity of the model has inspired the 
creation of a number of related models which have been amenable to 
rigorous combinatorial and probabilistic techniques. Most of the stud- 
ies, by far, have been on the probabilistic side, with systems being 
studied on the continuum, on infinite lattices in arbitrary dimension, 
or considering asymptotics of large finite lattices. The literature on 
Glauber dynamics in the Ising model is exceedingly large. When bonds 
between neighboring sites have random strengths, there seem to be a 
few studies for large or infinite systems, eg. [21 El lU IS]- We know of 
very few combinatorial results, eg [HI E] and they correspond to situa- 
tions where the bond strengths are not arbitrary. In what follows, we 
show that a model with Glauber-like dynamics exhibits rich algebraic 
and combinatorial structure. 

We consider a simplified one- dimensional version of the Ising model 
with spin-flip dynamics, on sites labelled 1 to L, and where there are 
two kinds of spins, which we label as occupation numbers {0, 1}. The 
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transition rules are given by an asymmetric version of the usual rule. 
Each site i looks only to the site to its left, i — 1, and the ith site 
switches its state with a rate that depends only on whether both sites 
are in the the same state or not. If the state of both sites is different, 
the transition occurs with rate oti (Glauber) and if they are the same, 
it occurs with rate (anti-Glauber). In other words, ajj's try to ho- 
mogenize the system and fy's try to make nearest neighbors opposite. 
In addition there are boundary interactions at the first site, which flips 
with rate ot\ if it contains a particle and f3\ if it does not. 

This model can be thought of as a voter model [8J within a hierarchi- 
cal demographic, that is to say individual % — 1 influences the opinion 
of individual i, but not vice versa. The rate at governs the tendency 
of i — 1 being a usual campaigner, whereas governs the tendency of 
i — 1 being a double agent. 

This model can also be thought of as a variant of the East model 
studied in [91 17] on the one dimensional lattice. In the east model, a 
spin was allowed to flip only if the site to its left was occupied and the 
rate of the flip depended on whether the site was occupied or not. 

Crisanti and Sompolinsky [3] considered a mean-field model of Glau- 
ber dynamics on Ising spins with asymmetric bonds with random stren- 
gths to understand properties of neural networks. This work is not 
related to it, but the idea is similar. In their work, asymmetry refers 
to partial asymmetry, instead of the total asymmetry assumed here. 

In an earlier work, the steady state of the border process for this 
model with = and ctj = 1 for i > 0, the asymmetric annihilation 
process, was obtained [10J by using a transfer matrix ansatz. The 
spectrum of the transition matrix for aij unequal was calculated in [TT] . 
This model, even with = is a clear generalization because one can 
distinguish here between particles and holes, whereas the former could 
not. In addition, there is no analog of the rate for in the earlier 
works because that would have amounted to creation of particles in 
the bulk. 

We will write down the precise statement of the Markov chain in 
Section |5J The main result of the paper is an explicit formula for 
the eigenvalues of the Markov chain for any size L, which we prove in 
Section [3J This leads, in an obvious way, to a formula for the spectral 
gap of the system. 

We will show that the transition matrix of the model exhibits a sim- 
ple recurrence relation allowing one to express the matrix for a system 
of size L in terms of that of size L — 1 in Section HI We prove a formula 
for the density in this disordered model in Section [5j We conjecture 
a formula for the normalization factor for the distribution, informally 
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called the "partition function" in Section El Lastly, we consider two 
special cases, which we naturally label the ferromagnetic and the anti- 
ferromagnetic limits in Section [7J 

2. The Model 

We consider a nonequilibrium system on a finite lattice with L sites 
labelled from 1 to L. Each sites is occupied by one of two spins 1 and 
0. The evolution rule in the bulk, for site i depending on the spin at 
site z — 1, is given by 

10—7-11 and 01 — > 00 with rate a,- , 

(1) 

y J 11 ^ 10 and 00 ->• 01 with rate ft . 

The evolution of the first site is given by 

1—^0 with rate ol\ , 

(2) 

y J 0^1 with rate ft. 

The rule at the left boundary fl2]) is constructed by supposing that 
there is a virtual site labelled to the left of the first site which is always 
empty. The left boundary conditions are deduced from the bulk rules 
(CO) by looking at the second component of the bond. 

Remark 1. There are two symmetries of the model; the first is mani- 
fested by interchanging all at 's and 's as well as flipping spins at all 
odd sites, and the second, by leaving a± and ft as they are, interchang- 
ing all other oti 's and 's and flipping spins at all even sites. 

Remark 2. This model has the property that correlation functions 
(tj^ . . . f]i n ) does not depend on the state of r/ in+1 simply because site 
i n + 1 cannot influence any site less than or equal to % n Further, as is 
well-known for Glauber dynamics, correlation functions of n sites like 
the one above depend only on those of sites less than or equal to n. 

3. Spectrum of the transition matrices 

We observe that the characteristic polynomial of the transition ma- 
trices factorize into linear factors and has an explicit formula. This 
is reminiscent of the formula for the asymmetric annihilation process 
conjectured in [10] and proved in [IT] . Let be the set of binary 
vectors of size L, that is vectors of length L whose elements are either 
or 1. 
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Theorem 1. The characteristic polynomial of Ml is given by 



(3) 



i=l 



For example, when L = 2, we have 
(4) 

|M 2 - Al 2 | = A (A + «i + fa) (A + a 2 + fa) (A + a x + a 2 + fa + fa) 

The proof of Theorem [I] will turn out to be a simpler version of the 
proof of the eigenvalues in the asymmetric annihilation process [TT] . 
We will first define a slight rearrangement of a Hadamard matrix. Let 
Bx be ordered lexicographically. For b, c € Ml, define the square matrix 
H L of size 2 L by 

(5) H L = ^ f(-l) Sc 



2 L/2 



where 6 is the reverse of b. Note that Hl is symmetric because 



(6) (Hl\ c = (-l) bc = (-l) b ' c = (H L ) c , b . 

This is different from the usual definition in which elements of Hl 
are written as (— l) 6 ' c . For example, the matrix for L = 3 is given by 

/ 1 
1 



(7) 



2V2 



1 


1 


1 


1 


1 


1 


1 \ 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


-1 / 



Lemma 2. H\ 



1/ 



Proof . We first look at the diagonal terms in H\. 
(8) {Hl\ b =Y,{H L \ c {HLU ' 

cgB L 



2 L 



\2b-c 



1. 



Now suppose b 7^ d. Then 



(9) 



L)b,d 



cen L 



cea L 
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where the addition of binary vectors is the usual xor-addition: x+ (1 — 
x) — 1 by definition and x + x = if x is a bit. 

Now, for any fixed b, d, we will construct an sign-reversing involution 
on Ml. Since b ^ d, b+d contains at least one entry equal to 1. Consider 
the leftmost such entry, at position i, say. To any element c, associate 
another element 



C f-> C = (Ci, . . . , C L -i, 1 - CL-i+X, CL-i+2, 



,c L ) 



Clearly, this is an involution and moreover, 



0. 



We have therefore partitioned Ml, for any fixed b,d,b ^ d, such that 
the number of terms which contribute +1 to the sum is exactly the 
same as that which contribute -1. Hence = 0. □ 

We now write down the transition matrix Ml in this language. We 
use the convention that (M^)^ represents the transition rate from 
c — > b, with the end result that M l \vl) = for the steady state column 
vector \vl). 



(10) (M L ) btl 



a, 



Pi, 



b + c = 0...0^0...0 

i 

and Cj_i + Ci = 1, 
b + c = 0...0^1^0...0 

i 

and Cj_i = Ci, 



E 



Ci; 



E ft. 



i=l 
Ci-i+Ci=X 



8=1 

i-l=Cj 



0, 



otherwise, 



Lemma 3. HlMlHl is lower-triangular. 



Proof. We will show the result by explicitly demonstrating the struc- 
ture of non-zero terms. By definition, 



[HlMlHl] 



c.d 



w E (-D 



c-e-\-f-d 



{M J 



L)e,f- 



eJ&L 



We will divide the proof into three parts. We first consider diagonal 
elements, then certain off-diagonal elements which are non-zero, and 
finally show that all other off-diagonal elements are zero by introducing 
another sign-reversing involution just as in Lemma [2j 
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C = d. 

Then c ■ e + f ■ d — c • (e + /). Suppose e differs from / only 
at position i, then (M L ) e j contributes either ctj of and the 
sign is always the same: +1 if q = and —1 otherwise. On the 
other hand, if e = f, this sum is identically and therefore we 
get a contribution of — ctj if + /j = 1 and — $ if /»_! = /j 
from (flDl) . 

To summarize, we get 2 i_1 terms contributing ±aj,±/9j, all 
of the same sign, from summands e ^ / and 2 i_1 terms con- 
tributing — ctj, — /3j from summands e = f. These will clearly 
cancel if q = and will give 2 L_1 (— 2a j — 2/3j) if q = 1. Thus, 

L 

(H L M L H L ) CiC = -J2ci{ai + A). 

»=i 

c^d 

For now, we look at conditions on c, d such that {Mi) e j = a, 
and (— l) c e+ / a! always has the same sign. This happens when e 
differs from / only at position % and + = 1 =>• ej_i = q. 
Now, 

c-e + f- d— e-c+f-d, 

= e ■ (c + d) + (e + /) • d, 
= /-(c + t?) + (e + /).e. 

Since e + / is fixed, only the former term contributes to the 
change of sign. As we vary e, this term does not change sign 
if and only if (c + d)j_i = (c + d)i and Cj = dj for all other 
values of j. Since we are assuming c ^ d, the only way this can 
happen is if c + d = ... 0110 ... with Is in the i — 1th and 
zth position. If this is the case, the sign of all these terms is 
determined just by dj, or equivalently, q. This gives us a total 
contribution of 2 L ~ 1 a i (— l) di . 

We also have contributions of on from diagonal terms e = / 
and q_i + q = 1 from ffTUj) . In that case, c • e + / • d = e.(c + d). 
Since we just argued that (c + d)j_i = (c + d)j, diagonal terms 
always contribute a term of (— 1) 1 (— a»), giving us a grand total 
of2 z - 1 a i . 

Therefore, diagonal contributions will cancel off-diagonal ones 
if di = 1, q = and will add to them if dj = 0, q = 1. 

By a very similar reasoning, we will get a contribution of — 
under exactly the same circumstances. The relation (c+d)i_x = 
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(c + d)i remains the same, but the roles are of c and d are 
reversed, giving us an overall negative sign. 
Assuming all other off-diagonal terms are zero, 

(H L M L H L ) c4 = on-Pi 

if and only if c and d differ only in the i — 1th and ith position 
and moreover di = and q = 1. Therefore d < c in the 
lexicographic order and the matrix is lower triangular. 

• Involution. 

To complete the proof, we have to show that when c and d 
do not satisfy c + d = . . . 0110 ... with Is in the i — 1th 
and ith position, (HlMlHl)^ = 0. To do this, we introduce 
another sign-reversing involution just as in Lemma [2J For a 
fixed c, d and any pair (e, /) which gives a contribution of a to 
(H L M L H L ) Cjd , say, we find another pair (e', /') which also give 

the same contribution but such that c-e + / -d — c- e' + f ■ d+1. 
Notice that (e, /) satisfies the condition that e and / differ only 
in the ith position and that /j_x + fi — 1. 

This is easily done. Since c ^ d, find the smallest value of j, 
j 7^ i, i — 1 such that Cj + dj = 1. Then choose (e', /') as 

e H*e' = (d, . . . , e L -j, 1 - e L - j+1 , e L _ j+2 , . . . , e L ), 

f = • • • ) /i-j) 1 — fh-j+l, Jl-j+2, • • • j /l)- 

Then e + f = e' + f and e.(c + d) = e'.(c + d) ± 1 and we are 
done. Notice that we need j ^ i — 1, i to ensure that ej_i = e« 
remains intact. If c differs from d only at these two positions, 
then we already proved that there can be no such involution. 

□ 

We have now done all the work necessary to write down a proof of 
the main result. 

Proof of Theorem [Tl Since Hl is symmetric and H\ = 1 by 
Lemma H^M^H^ has the same characteristic polynomial as M^. 
But the former is lower triangular. The characteristic polynomial is 
therefore 

(15) \M L - A1 L | = n(" A " (H L M L H L ) iti ). 

i=l 

Using the representation of the rows and columns using boolean vectors 
of size L and using (1121) . we have the required result. ■ 
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Corollary 4. The spectral gap for the system of size L is given by 



L 



(16) min(aj + 



i=l 



Proof. The proof directly follows from the explicit construction of eigen- 
values in Theorem [TJ □ 

4. Recurrence Relation for the Transition Matrices 

We will prove that there is a recursion of order one among the tran- 
sition matrices. Let M L (ai, . . . , «£,; fti, . . . , (3l) represent the transition 
matrix for the system of size L with parameters ctj and 

Theorem 5. Let 1 denote the identity matrix of size 2 L ~ l . Then 
can be expressed in 2 x 2 block- diagonal form as 

(17) M L {a u ...,a L ',p x ,...,p L )-- 

where 



M 1A 


ail 




M 2i2 



(18) 



M^i = M L _i(a 2 , «3, • • • , "l; #2, Ai, • • • , - Al, 
M 2 , 2 = M L _i(^2,«3, a 2 ,^3, ail, 

with the initial matrix for L = 1 pnien by 

(19) M 1 (a 1 ;A)=(^ 1 

Proof. The matrix for size 1 given in ( !T9|) is easy to check. The matrix 
recursion ( flTl) is proved by looking at transitions that binary vectors 
beginning with different two-bit vectors undergo. Let v, v' denote bi- 
nary vectors of length L and w, w' denote vectors of length L — l. There 
are then four possibilities of transitions depending on the first bits of 
v and v', namely v = Ow or v = lw and similarly for v' . 

(1) Ow — » lw' (resp. lw —> Ow'): 

These transitions are only possible if w = w' and then oc- 
curs with rate fli (resp. «i). This is why the (2,1)— th (resp. 
(l,2)-th) block of M L is ftl (resp. a x l). 

(2) Ow Ow' (resp. lw —> lw'): 

In this case, all the transitions for the system of size L — 1 go 
through but with rates whose indices are one more than that 
of the smaller system simply because one extra site has been 
added to the left. For example, what occurred with rate a 2 in 
the system of size L — l occurs with rate now. Because of the 
extra transitions in the off-diagonal blocks (described above), 
we must add negative terms to the diagonal so that column 
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sums are zero. Therefore we subtract f3\l (resp. ail) from 
M M (resp. M 2>2 ). 

□ 

For later consideration, we also express the recursion relation using 
a second order recurrence relation. Expressing the transition matrix 
for the system of size L — 1 as 

M L (a 1 ,...,a L ;p 1 ,...,^ L ) = 
( Mn- a 2 l ail \ 



(20) 



(Pl + & 


)1 


a 2 l 


ail 





/3 2 1 


M 2j2 - 
(A + a 2 )l 





ail 


fill 





Mi,i- 
(ai + a 2 )l 









/?il 


a 2 l 


M 2>2 - 
(«i + /3 2 )l 


Mi,x = 




_ 2 (a 3 , a 4 , . . 


, Oil] /?3, 04, ■ 


■■,&), 


M 2i2 = 




-2(^3, a 4 , . . . 


, «x; a 3 , /3 4 , . 


■-,&)■ 



where 
(21) 



Remark 3. Notice that the transition matrix Mi is invariant under 
the transformation M L -B- (M L ) T and a •<->• /3. 

5. Density in the Steady State 

We compute the density of l's in the steady state at site k in a 
system of size larger than k. Using Remark [21 this density is completely 
independent of the size of the system. 

Let Sk be the set of subsets of [k] — {1, . . . , k} with an odd number 
of elements, and for s e Sk, let s = Sk \ s. 

Lemma 6. The density of 1 's at site k is given by 

(22) (m) 



sGSfe i&s j£s 



Ufa + to 

3=1 

Proof. The master equation for the density of l's at site 1 is 
d 



(23) 

which implies 771 



dt 



(Vi) = - Vi) ~ txiim) = °> 



(ai+ft)' 



which is consistent with the 



with k — 1. 
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For k > 1, the master equation is 

( 24 ) ^(Vk) = - %-i)(l - %)) + a fc (r/ fc _i(l - 7/ fc )> 

- a fc ((l - r/ fc _i)?7 fe ) - (3 k (Vk-iVk) = 0, 
which on simplifying gives 

(25) (77fc) = ——-5 . 

«fe + Pk 

One can check that 0221) satisfies this recurrence. □ 



6. The Normalization Factor 

Consider the system of size L. For any configuration 77, we can 
express its steady state probability as a numerator divided by a de- 
nominator. We define the normalization factor Zl{o.,P) as the least 
common multiple of the denominators of the steady state probabilities 
for all configurations n such that the greatest common divisor of the 
corresponding numerators is 1. 

We have been able to find a conjecture for Zl but have not been 
able to prove the formula in general. In certain special cases discussed 
in Section [7J the following conjecture does reduce correctly. 

Conjecture 1. The normalization factor for a system of size L with 
rates cti and is given by 

L 

(26) Z L (aJ)=H(a i + p i ) J] + & + % + 

i=l l<i<j<L 



7. Special Cases: Ferromagnetic and Antiferromagnetic 

MODELS 

We now analyze two limits of this model which depend on two pa- 
rameters, a.\ = a and /3i = /3. Both these limits can be interpreted 
as nonequilibrium analogs of the ferromagnetic and antiferromagnetic 
Ising models respectively, hence the terminology. 

The strategy for the analysis of both cases is similar. We will first use 
Theorem [S] to write a first order recurrence for the transition matrices. 
We will then show that the transfer matrix ansatz [10] holds, using 
which we will express the steady state probabilities of the system of 
size L in terms of the system of size L — 1. 

We first recall the definition of the transfer matrix ansatz. A family 
of Markov processes satisfies the transfer matrix ansatz if there exist 
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matrices Tl for all sizes L such that 

(27) M L+l T L = T L M L . 

We also impose that this equality is nontrivial in the sense that 

(28) M L+1 T L ^ 0. 

The rectangular transfer matrices or conjugation matrices Tl can be 
interpreted as a consequence of the integrability of the model. 

For the system we consider here, the transition matrices are of size 
2 L . It is most convenient for us to take the naturally ordered basis of 
binary sequences of size L. For example, when L = 2, the ordered list 
is (00,01,10,11). 

Using the transfer matrices, we will calculate the normalization fac- 
tor Zl for both these models. Letting \vi) be the Perron- Frobenius 
eigenvector of the system of size one, we can define 

(29) Z L = (l L \T L -iT L . 2 .--T 1 \v 1 ), 

where (1l\ is the row vector of size 2 L composed of l's. 

For these special cases, we will explicitly calculate the density, the 
recurrences for the transition matrix and the transfer matrix, as well 
as the formula for the normalization factor. We will omit all the proofs 
since they are quite straightforward, if somewhat tedious. The recur- 
rences for the transition matrix follow from the recurrence (120 p and 
the others can be provided without too much difficulty using induc- 
tion. The technique of proof is the same as was used in [TO*]. 



7.1. Ferromagnetic limit. We consider the case ot{ = 1 and = 

for i > 1. The densities are easily calculated. The only term that 
contributes to the sum in ( T22|) is the subset {1} G Sk- Therefore, 

(30) (%) 



for all k. Therefore all densities are identical, which is reminiscent of 
the ferromagnetic limit where all spins prefer to be aligned in the same 
direction. 

Corollary 7. Let t denote the identity matrix of size 2 L ~ 2 . Suppose 
the transition matrix for size L — 1 is written in block- diagonal form as 

(31) M[ F \ - 





al 


01 
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Then 



(32) 



M 



(F) 



( m[? 


i 


al 








(l-a + P)l 





al 







Mjf- 
(1 -£ + c*)l 





\ 


f3l 


1 


Mi? 



/ 



where M L is written as a 2x2 fr/ocA; matrix with each block made up 
of matrices of size 2 L_1 . The initial matrix for L — 1 is given by 



(33) 



M[ F) = 



-/3 a 
(3 —a 



Let gl denote the matrix of size 2 L with l's on the antidiagonal and 
zeros everywhere else. For example, 



(34) 



o-i 



r 

1 or 



Lemma 8. Define rectangular matrices T^ of size 

2 l+i 

x 2 L using 

the following recurrence relation. 



(35) 



7f> = 



/ 1 — /3 + a a \ 

a 

f3 

\ /3 l-a + / 



and define T^ by the following recursion. If 



(36) 



T 



(F) 



L-l 



/ T {F) 
1 1 11 


1 12 


\ 





T (F) 
1 22 


T (F) 
J 31 







J 42 


/ 
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then 



(37) 



T, 



( wi? 


T (F) 
1 11 


9T (F) 
ZJ 12 


1 12 \ 










T (F) 
1 22 








T (F) 
1 12 











T {F) 


0T (F) 
ZJ 22 


9T (F) 


T (F) 
J 31 











T (V) 
1 Al 








T (F) 





rp(F) 

°~L-2-l-42 °~L-2 







rp{F) 
ZJ 41 


T iF) 
1 i2 


2Tif J 



Then the family of Markov processes given by the transition matrices 



{M^ } satisfy the transfer matrix ansatz with T£ 



(F) 



defined above, ie 



L-l 



T (F) M (F) 



Using Lemma[8l we can compute the normalization factor Z^\a, j3) 
using fl29l) by setting 

(38) \ Vl ) = 

For example, Z\ = a + f3. We find a remarkable property of the 
partition function of the system, namely its super-extensive growth 
with the size of the system. 

Corollary 9. The partition function of the system of size L is given 
by 



(39) 



Zf ] = 2(^ 1 )(a + /3)(l + a + /?) 



L-l 



7.2. Antiferromagnetic limit. We consider the case «j = and /3j 
1 for i > 1. Now the only terms that can contribute to the sum in 
are the subset s — <f), {1} G S^. Therefore, the density depends on the 
parity of the site. If k is odd (resp. even), s = [k] (resp. s — {2, . . . , A;}) 
is the only contribution, 



(40) 



(Vk) 



I - 9L - 



k odd, 
k even. 



This is reminiscent of the antiferromagnetic Ising model in which neigh- 
boring sites prefer to be aligned opposite to one another. 
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Corollary 10. Let 1 denote the identity matrix of size 2 L 2 . Suppose 
the transition matrix for size L— 1 is written in block- diagonal form as 



(41) 





al 


/31 





Then 



(42) 



M 



(A) 



( M[f> - 1 





al 


\ 


1 


m£>- 

(a-/3)l 





al 


PI 





M {A) - 
(f3 - a)l 


1 


\ 


/31 





m£> - 1 J 



where is written as a 2 x 2 block matrix with each block made 

of matrices of size 2 L ~ 1 . The initial matrix for L = 1 is given by 



up 



(43) 



(A) 



—(3 a 
(3 -a 



Lemma 11. Define rectangular matrices T^ of size 

2 l+i 

x 2 L using 

the following recurrence relation. 



(44) 



T (A) = 



( a 

1 — (3 + a a 
(3 1-a + p 
13 



and if 



(45) 



T, 



(A) 



L-l 



( o 



41 



T (A) \ 
1 12 » 



T (A) 


rp(A) 


rp(A) 

J 31 


rp(A) 



o / 
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then 



(46) 



T, 



(A) 












1 A'\ 
rp(A) 

1 12 








0T {A) 
z ' 1 22 


rp(A) 
1 22 




2aT^fa 


rp(A) 
1 12 


r\rp{A) 

Z ' 1 12 


T (A) 

1 2l 





T (A) 

1 22 








T {A) 
J 31 





rp[A) 

1 32 


cyrp(A) 

ZJ 41 


T {A) 






T (A) 

J 31 


cyrp(A) 
Z1 31 








T (A) 








o 



/ 

Then the family of Markov processes given by the transition matrices 

^ defined above, ie 



{M^ } satisfy the transfer matrix ansatz with 



MjT, 



L-l 



rp(A) M (A) 

L L-\ m L-\- 



Using Lemma[Hl we can compute the normalization factor Z^''(a, /?) 
using fl29l) by using the same formula for \vi) as in the ferromagnetic 
case, (1381) . We find that the normalization factor is exactly the same 
as for the ferromagnetic model. 

Corollary 12. The partition function of the system of size L is given 
by 



(47) 



r(A) 



2^){a + (5){l + a + (5) L -\ 



7.3. Relation between these special cases. The transition matrix 
M L for X = F, A is invariant, in addition to the transposition sym- 
metry, under the transformation 

(48) (M { L X) ) i;j O (Aff V+i-i^+w and a O /? 

in the usual matrix coordinate notation where i,j vary from 1 to 2 L . 
Similarly, the transfer matrix T^ is invariant under the transforma- 
tion 

(49) (Tf \ 3 <-> (Tf ) ) 2 , +1+1 _ li2 , +1 _, and a <-> (5 

where this time i varies from 1 to 2 L+1 and j varies from 1 to 2 L . The 
normalization factors for the two systems are exactly the same, Z L = 
. One also sees that the recurrence relations for the transition 
matrices and the transfer matrices are similar. 

The reason that there are so many similarities between the two mod- 
els is that one can establish an exact correspondence between the fer- 
romagnetic and antiferromagnetic limits defined previously using Re- 
mark [TJ In this case, we use the second part of the remark. As a 
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consequence of this coupling, time-dependent correlation functions in 
the ferromagnetic model are related to those in the antiferromagnetic 
model. 
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